function model=update_sigma_01(model)

% global model;
% 

N         = model.N;

temp      = [sqrt((1/N)*sum(model.sigma_1.^2+(repmat(model.alpha_01,N,1)-model.alpha_1).^2,1))];
ind       = find(temp<model.MINVALUE);
temp(ind) = model.MINVALUE;

model.sigma_01 = temp; 

end
